Phase diagram of quasi-two-dimensional bosons in laser speckle potential 
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We have studied the phase diagram of a quasi-two-dimensional interacting Bose gas at zero tem- 
perature in the presence of random potential created by laser speckles. The superfluid fraction and 
the fraction of particles with zero momentum are obtained within the mean-field Gross-Pitaevskii 
theory and in diffusion Monte Carlo simulations. We find a smooth crossover from the superfluid 
to the insulating state, when the strength of the disorder grows. Estimations of the crossover lines 
are compared with the predictions of the percolation theory in the Thomas-Fermi approximation. 
Analytical expressions for zero-momentum fraction and superfluid fraction are derived in the limit 
of weak disorder and weak interactions within the framework of the Bogoliubov theory. Limits of 
validity of various approximations are discussed. 
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Introduction. — The interplay of disorder and interac- 
tions in quantum systems is one of the fundamental prob- 
lems in modern physics. Whereas disorder tends to lo- 
calize non-interacting particles, giving rise to Anderson 
localization and suppression of spatial coherence, weak 
repulsive interactions tend to bring the system back to 
the extended state without necessarily restoring the co- 
herence due to the enhancement of quantum fluctuations. 

Ultracold atoms provide an irreplaceable tool to study 
interacting quantum systems with disorder due to the 
unique opportunity to control the system [1]. Interaction 
strength can be tuned by Feshbach resonances, while dis- 
order with known statistical properties and tunable pa- 
rameters can be created either by optical means using 
incommensurate optical lattices and laser speckles [2] or 
through the interaction with impurities localized at ran- 
dom positions [3] . Using confining potentials in different 
spatial directions, one can also control the system dimen- 
sionality. 

Laser speckles allow to create unbounded random po- 
tentials with a finite correlation length [4, 5]. In recent 
years, laser speckles were used for the experimental ob- 
servation of Anderson localization of bosons in one di- 
mension (ID) [6] and three dimensions (3D) [7] as well 
as of fermions in 3D [8] , suppression of transport in elon- 
gated ID geometry [5, 9-11] and Bose-Einstein conden- 
sate in 3D optical lattices [12]. Coherent phenomena and 
diffusion were also experimentally studied in quasi-2D ge- 
ometry [13-16]. 

The properties of single-particle eigenstates in the pres- 
ence of laser speckles are almost the same as in the case of 
uncorrelated disorder. The difference was found only in 
ID, where there is a crossover from exponential localiza- 
tion to the algebraic one [6, 17], while all the eigenstates 
in the presence of uncorrelated disorder are exponentially 
localized. In 2D, the wave functions always show an ex- 



ponential long-range decay and one has to distinguish be- 
tween weak and strong localization [18, 19]. In 3D, there 
is a critical energy that separates extended and localized 
states (mobility edge) but its precise determination still 
remains a difficult problem [20]. 

The quantum many-body problem of interacting parti- 
cles in the presence of disorder is even more challenging. 
Most of the studies of the interacting bosons in the pres- 
ence of laser speckles are done within the framework of 
the mean-field theory. Gross-Pitaevskii equation (GPE) 
was used to study ground-state density profiles in ID [21], 
expansion in ID [17, 22-25], transport of quasi-lD BEC 
in the subsonic and supersonic regimes [26], dipole and 
quadrupole collective oscillations in 3D [22]. Thomas- 
Fermi approximation was successfully employed in the 
regime when the correlation length of the disorder po- 
tential is much larger than the healing length of the con- 
densate [5, 21, 24, 27, 28]. Combined with the percola- 
tion theory [27, 29], this approximation was used to esti- 
mate the critical values of parameters for the superfluid- 
insulator transition in a 3D-system [27]. 

Quantum fluctuations were taken into account within 
the framework of the Bogoliubov theory that predicts 
Anderson localization of quasiparticles [30, 31]. The Bo- 
goliubov theory was also used to study excitation spec- 
trum, corrections to the speed of sound and condensate 
population in the case of weak disorder [32-34]. GPE 
and Bogoliubov-de Gennes equations in ID reveal the 
supcrfluid-Bose-glass transition [35, 36], although in an 
earlier paper [37] it was argued that this is rather a 
crossover. Quantum Monte Carlo (QMC) calculations 
in 3D were performed in Refs. [27, 38], where the finite 
temperature superfluid-insulator transition was studied. 

In the present work, we investigate the superfluid- 
ity and the off-diagonal long-range order of interacting 
bosons at zero temperature in the presence of laser speck- 
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Figure 1. (color online) Phase diagram. Lines of crossover 
from the superfiuid to the insulator regime calculated from 
the requirement that the superfiuid fraction takes fixed values 
v s = O.f (blue), 0.5 (red), 0.9 (green). Dashed line is the 
estimate of the critical values of V$ from the percolation model 
in the local-density approximation [Eq. (12)]. The insets show 
density profiles calculated for one disorder realization with 
Vs/fkJho = 0.1, a s /a ho = 0.01 (a); V 3 /fiu)ha = 1, a s /aho = 
0.01 (b); V s /ftw bo = 1, as/aho = 1 (c); V s /fajJho = 40, a s /a ho = 
1 (d). 



les in a quasi-2D setup. On the basis of QMC calcula- 
tions we come to the conclusion that there is a smooth 
crossover from the superfiuid to the insulator phase and 
the crossover lines that determine the dependence of the 
disorder strength V£ on the interaction strength a s at 
fixed values of the superfiuid fraction v s are shown in 
Fig. 1. 

The system. — The Hamiltonian of N interacting 
bosons in the presence of an external random potential 
V(r) and a tight harmonic trapping in z direction is 



H 



N 

E 

i=l 



2m 



V(Ti) 



We assume that the frequency Who is so large, that parti- 
cles do not populate higher levels of the trap. Under this 
realistic assumption, the system is effectively quasi-two- 
dimensional. We consider a system of a finite size L x L 
in the (x, y)-plane. The interaction potential V pp (r) is 
chosen to be of a short-range type as a relevant model 
for ultracold atoms. In the QMC simulations, we model 
those as hard spheres of the diameter a s . In the dilute 
regime, the potential can be replaced by a contact term 
4irh 2 a s /m S(r) that is used in the GPE. For both inter- 
action potentials, a s is the s-wave scattering length. 

The external random potential V(r) is assumed to be 
created by laser speckles [4, 5]. In the presence of a strong 
harmonic confinement one can neglect its ^-dependence. 
We numerically generate two-dimensional speckle pat- 
terns V(r) = V(r±) following the procedure described 
in Ref. [39]. Each realization of V(r±) is a periodic func- 
tion of x and y, therefore, permits us to impose periodic 
boundary conditions in the (x, y) -plane. V(r±) has ex- 



ponential probability distribution with the mean value 
V s = V(rjJ that is an important parameter of the prob- 
lem. Explicit form of the autocorrelation function 



f(r± - r' ± ) = V(r ± )V(r' ± ) - V(r ± ) 



(1) 



is determined by the aperture of the laser system [4] . In 
the case of a square aperture and in the limit L — > oo, it 
is given by 



/(rj 



c' x ) = V s 2 smc 2 



y-y 

Lc 



(2) 



where L c is the correlation length that determines the 
characteristic size of the speckles. In the case of a circular 
aperture, it takes the form 



f(r^-r' ± ) = V s 2 



Ji(x \r± - r' ± \/L c ) 



(3) 



where xq ~ 3.83171 is the first nontrivial zero of the 
Bessel function J\{x). 

Numerical methods. — For the numerical solution of 
the problem we employ variational and diffusion QMC 
methods [40]. In the variational QMC calculations, the 
trial many-body wave function is constructed in the form 



N 



T (r 1 ,...,r N )=l[<S>(r l )Y[f 2 (\r i 



i—l i<j 

where /a(|r|) is a two-body Jastrow term. The strong 
harmonic confinement in the z direction allows to repre- 
sent the one-body term as $(r) = tp(r±)tpi 10 (z), where 
V'ho('Z) is the ground-state wave function of the harmonic 
oscillator and ip(r±) satisfies the two-dimensional GPE 



2m 



V(v A _)+g 2 N\^(v 1 _)Y 



(4) 



with the effective coupling constant g 2 = 
\f%Tih?a s / (maho)- 

The superfiuid fraction v s is obtained in QMC sim- 
ulations using the winding number technique [41]. It 
can be also obtained from the GPE in a moving refer- 
ence frame introduced by means of substitution ip(r±) = 
4>{r±) exp(ik s • rj_) in Eq. (4), where </>(rj_) satisfies peri- 
odic boundary conditions. For small k s , the expectation 
value of the momentum P has a linear dependence on 
k s , i.e., P = Nv s fik s which determines the superfiuid 
fraction v s . 

Fraction of particles with zero momentum isq which 
gives a lower estimate of the condensate fraction [42] is 
calculated from the one-body density matrix (OBDM) 
/?i(r,r'), integrated over z' — z, as the asymptotic value 
at large distance |rj_ — TjJ — > oo. In QMC calculations, 
OBDM is obtained by extrapolation from the variational 
and mixed estimators, while in the mean-field theory, 
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Figure 2. (color online) Dependence of the superfluid fraction 
on the system size for a s /aho = 0.1. Solid symbols indicate 
the results obtained from the GPE; open symbols with DMC; 
Straight lines are the fits (5) corresponding to (from top to the 
bottom) V s /huj ho = 0.1, 0.18, 0.32, 0.56, 1, 1.8, 3.2, 5.6, 10. 
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Figure 3. (color online) Superfluid fraction in the thermody- 
namic limit. The results obtained from the GPE are shown 
by solid lines, QMC data are indicated by symbols, and the 
predictions of the Bogoliubov theory described by Eq. (8) are 
represented by dashed lines. Different colors distinguish be- 
tween the values of a s . 



it is readily obtained from the solution of the GPE as 
Pi(r,r') = *(r)$(r'). 

Numerical results. — In the problem we consider, there 
are three natural spatial scales: healing length £ = 
fi/^/rng 2 'T-2, correlation length of the disorder L c , and 
the oscillator length a\ lo = y/H/ (mwho)- We use ctho as a 
length unit and the energies are measured in the units of 
hujh . Our numerical calculations are done for the square 
aperture with realistic values of n 2 = N/L 2 = 0.1/a^ o 
and L c = ai lo . 

The density profiles calculated for one disorder realiza- 
tion in different regimes are shown in Fig. 1. In the case 
of weak interactions [insets (a) and (b)], corresponding to 
£ rs 14 Oho larger than L c , a rather weak disorder is suf- 
ficient to localize the atoms in the region where the ran- 
dom potential has an extended minimum. In the case of 
strong interactions [insets (c) and (d)], where £ rs 1.4 Oho 
is comparable to L c , the density profile has a multipeak 
structure with maxima corresponding to local minima of 
the disorder potential. We find that the density profiles 
are very well described by the GPE even for very strong 
disorder similarly to the 3D case [27]. 

The superfluid fraction calculated in finite-size systems 
for a s /aho = 0.1 and different values of V s after averaging 
over several disorder realizations is shown in Fig. 2. It 
decreases with the increase of V s and L. The size depen- 
dence (see Fig. 2) is well approximated by the fit 



v s (L) = v s {oo) + c/L 



(5) 



where 14(00) is the value of the superfluid fraction in 
the thermodynamic limit and c describes finite-size cor- 
rection. The dependence of ^ s (oo) on is reported in 
Fig. 3. The predictions of the GPE are in good quantita- 
tive agreement with the QMC data as long as the inter- 
actions are not too strong. Quantum fluctuations fully 
taken into account in the QMC calculations and com- 



pletely neglected in the GPE become noticeable for large 
interactions and lead to additional suppression of the su- 
perfluidity. Our calculations show that the superfluid 
fraction becomes larger for stronger interactions. We also 
observe smooth decrease of v a to zero when the disor- 
der strength grows. This leads us to the conclusion that 
there is a crossover from the superfluid to the insulat- 
ing regime rather than a phase transition. The crossover 
lines V s = V s c (et s ) for fixed v s are shown in Fig. 1. 

The fraction of particles with zero momentum shows 
almost no size dependence for large enough L. The ther- 
modynamic values ^o(oo) presented in Fig. 4 are obtained 
by fitting the GPE data for N > 40 by a constant. With 
the increase of disorder, vq smoothly decreases similarly 
to iv For stronger interactions, the GPE predicts an 
increase of vq, while exact QMC calculations show in 
general a non-monotonic dependence of vq on a s . The 
comparison of the results obtained from the GPE and 
QMC shows that quantum fluctuations play an impor- 
tant role for a s > 0.5 Oho and lead to the suppression of 
v Q . 

Bogoliubov theory. — For weak external potential V(r), 
the GPE equation can be solved analytically in the lowest 
orders of V [21, 31-34, 42]. Assuming that V(r) has a 
period L in d spatial dimensions and possesses cubic or 
spherical symmetry, perturbative solution of the GPE in 
the moving reference frame in the limit of vanishing k s 
leads to the following result for the superfluid fraction: 
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where gd is an effective interaction parameter in d dimen- 
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Figure 4. (color online) Fraction of particles with zero mo- 
mentum in the thermodynamic limit. The results obtained 
from the GPE are shown by solid lines, QMC data are in- 
dicated by symbols, and the predictions of the Bogoliubov 
theory are represented by dashed lines. Different colors dis- 
tinguish between the values of a B . 



sions, rid = N/L d , 
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Within the same formalism we also get vq = 1 — j ( 1 — v s ). 
Standard Bogoliubov theory [43] gives exactly the same 
result for v s but Vq contains also an additional depletion 
—g2m/(4nh 2 ) [33, 34] due to quantum fluctuations [44] 
which are neglected in the GPE. 

Since the speckle potential is random, we have to per- 



form statistical averaging in Eq. (6). In general, V(n 

coincides with the Fourier transform (7) of the autocorre- 
lation function (1) of the disorder potential. In the case 
of (5-correlated disorder, this allows us to reproduce the 
results of Ref. [43]. For laser speckles with the square 
aperture and in the limit L — > oo, we get 
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(8) 



where b = L c /(7r£). For the circular aperture, we have 



v* = 1 - 2 



Vsb 

52 "2 



2b 2 - 2b\/l + b 2 + 1 



(9) 



where b = L c /{xq^). In the limit £ 3> L Cl v s tends to 
unity, while in the opposite limit it takes the asymptotic 
value v s = 1 - i (V s /g2n 2 ) 2 . 

The predictions of the Bogoliubov theory for v s and 
Vo are shown by dashed lines in Figs. 3, 4, respectively. 
As expected Bogoliubov theory works better for stronger 



interactions, although it has a tendency to underesti- 
mate v s as well as v§ already at small disorder strengths. 
Higher order calculation would probably improve the sit- 
uation. 

Thomas-Fermi regime and percolation analysis. — If 
£ <C L Cl the density profile of the system in an exter- 
nal potential can be calculated within the Thomas-Fermi 
approximation and has the following form 



V(r x ) 



g 2 N 
g 2 n 2 + V s 



e[/z-F(rx)], (10) 
-V B exp(-fi/V a ) . (11) 



Equation (10) allows the existence of the regions where 
the density completely vanishes. The fraction of space 
occupied by the particles is O = 1 — exp(— fi/V s ). If Q 
exceeds the critical value of the percolation transition J7 C , 
there are infinitely extended regions of space occupied by 
the particles and the system is superfluid. Otherwise, it 
is insulating. Such reasoning leads to the conclusion that 
there must be a sharp transition from the superfluid to 
the insulator which is in contradiction to the results re- 
ported above. Nevertheless the percolation analysis can 
be used for the estimation of the "critical" values of pa- 
rameters as this was done in 3D where there is indeed a 
sharp transition [27]. Using Eq. (11) we get an estimate 
for the critical value of the disorder strength 



K c = - 



52 n 2 



n c + ln(l - fi c ) 



(12) 



Numerical study for two-dimensional systems gives tt c « 
0.407 [27, 29]. This result is shown in Fig. 1 by the 
dashed line and predicts qualitatively correct behavior 
of the crossover lines. 

Conclusions. — Using mean-field theory and QMC cal- 
culations we have studied the superfluid fraction v s as 
well as the fraction of particles with vanishing momen- 
tum pq for the interacting Bose gas at zero temperature 
in the presence of correlated disorder created by laser 
speckles in a quasi-2D setup. It is shown that v s and 
vq smoothly decrease with the disorder strength V s indi- 
cating a crossover from the superfluid to the insulating 
regime. With the increase of interactions a s , v s always 
grows while vq has in general a non-monotonic depen- 
dence. For large a s and V s , vq can exceed v s . The results 
of the mean-field theory are in excellent agreement with 
the exact QMC calculations in the dilute regime even in 
the case of strong disorder. The predicted phase diagram 
can be accessed in the present-day experiments. 
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